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A  very  important  sub-problem  in  the  task  assignment  problem  for  unmanned  air  vehicles 
(UAVs)  is  the  evaluation  of  costs  for  the  state  transitions  of  a  directed  graph.  Usually 
a  Dubins  vehicle  flying  in  the  absence  of  wind  is  considered  in  the  computation  of  costs. 
However,  when  a  prevailing  wind  vector  field  is  considered,  the  costs  taken  on  very  different 
values  and  the  task  assignment  problem  can  have  very  different  solutions.  In  this  paper, 
we  consider  the  problem  of  constructing  minimum  time  trajectories  for  a  Dubins  vehicle  in 
the  presence  of  a  time  varying  wind  vector  field.  We  present  results  on  the  existence  and 
uniqueness  of  minimum-time  solutions  for  a  Dubins  vehicle  flying  in  a  general  time-varying 
wind  vector  field  under  some  technical  conditions.  These  results  extend  the  conclusions  of 
the  well-known  Dubins  theorem.  We  also  propose  an  algorithm  for  obtaining  the  minimum¬ 
time  solution  for  an  UAV  and  prove  its  convergence.  We  also  present  the  results  of  numerical 
experiments  that  show  that  the  importance  of  considering  wind  vector  fields  while  planning 
the  tour  for  an  UAV. 
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Nomenclature 

Unmanned  Air  Vehicle 

Micro  Air  Vehicle,  usually  with  wing  span  approximately  1  ft. 
Constant  Speed  of  the  MAV. 

Time  in  seconds. 

Free  final  time  for  the  trajectory  planning  problem. 

Position  of  the  center  of  mass  of  the  MAV  in  IR2. 

Heading  of  the  MAV  (in  radians). 

=  Ou# 

Initial  position  and  heading  at  initial  time  0. 

Desired  position  and  heading  at  final  time  tf. 

Steering  control  input  to  the  MAV. 

Upper  bound  on  \u(t)\  for  all  t. 

The  wind  vector  field  is  a  function  of  x  and  t  in  general. 

Angle  by  which  the  initial  coordinate  system  is  rotated  for  numerical 
stability  of  the  trajectory  planning  algorithm. 

Position  of  the  MAV  after  a  coordinate  transformation  for  numerical 
stability. 

Heading  of  the  MAV  after  a  coordinate  transformation. 
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Introduction 


Cooperative  Control  of  multiple,  autonomous  unmanned  air  vehicles  (UAVs)  is  an  active 
area  of  research  that  holds  enormous  potential  for  military  and  civilian  applications.1'3  This 
new  paradigm  for  control  has  been  implemented  in  the  MultiUAV  simulation  platform  by 
the  Air  Force  Research  Laboratory.3  It  has  a  hierarchical  architecture,  where  at  the  highest 
level  the  dynamics  of  the  controlled  agents  are  usually  suppressed  and  a  task  allocation  for 
the  agents  is  performed  using  graph  theory.  The  tasks  for  the  agents  are  usually  tightly 
coupled  in  time3,4  and  hence  estimates  of  the  times  taken  for  each  agent  to  fly  from  one 
destination  to  the  next  is  highly  critical  for  correct  assignment  of  tasks.  This  estimate  of 
times  is  usually  obtained  by  considering  a  kinematic  model  of  the  air  vehicle,  along  with 
kinematic  turn-radius  constraints,  to  keep  the  computation  time  at  a  manageable  level.5'' 
The  key  result  that  is  used  in  this  computation  is  Dubins’  result  on  the  existence  of  minimum 
time  solutions  for  a  kinematic  model  with  minimum  turn-radius  constraint.8  However,  this 
result  is  only  for  valid  for  zero- wind  and  hence  all  of  the  available  cost  estimation  algorithms 
are  only  valid  for  zero-wind.  In  this  paper,  we  propose  a  new  method  for  the  computation  of 
minimum  time  solutions  for  the  trajectory  planning  problem  for  a  Micro  Air  Vehicle  (MAV) 
or  UAVs  in  the  presence  of  time-varying  wind,  and  prove  its  convergence.  We  also  present 
the  results  of  numerical  experiments  that  show  that  the  order  in  which  targets  need  to  be 
over-flown  to  minimize  total  flight  time  can  be  very  different  if  the  winds  are  accounted  for. 

We  model  an  MAV  flying  with  a  constant  speed  in  the  wind  axes  and  at  a  constant 
altitude.  The  kinematic  equations  of  motion  for  the  MAV  are: 

x  =  U(cos,0,sin-?/>)  +  W(x,t)-,  ip  =  u.  (1) 

where  x  =  (x^x-p),  W(x,t)  =  {W\{x,  t),  W2(x,  t)),  and  Wi(x,t)‘,  i  —  1,2  are  functions  with 
bounded  derivatives.  These  equations  contain  the  wind  vector  field  that  is  not  considered  in 
earlier  works.5'9  Of  course,  the  actual  model  for  a  MAV  consists  of  Euler’s  equations10  and 
kinematic  model  used  above  implies  a  perfect  response  to  the  turn  commands.  However,  it 
must  be  kept  in  mind  that  the  kinematic  model  is  only  employed  to  compute  the  order  in 
which  the  targets  are  to  be  visited.  Once  the  order  is  determined,  optimal  control  theory 
should  be  employed  to  determine  the  actual  time-optimal  path  to  be  flown  by  the  aircraft. 
A  major  reason  for  using  the  kinematic  model  is  the  fact  that  only  necessary  conditions  for 
optimality  exist  for  the  second  order  model  (given  by  Pontryagin’s  Maximum  Principle)  and 
no  results  on  existence  of  solutions  are  available.  The  situation  is  different  for  the  kinematic 
model  where  in  the  absence  of  wind  (that  is  W(x,t)  =  (0,0)),  the  well-known  Dubins’ 
theorem8  posits  the  existence  of  a  time-optimal  solution  for  any  combination  of  initial  and 
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final  positions  and  velocities  of  the  UAV  or  MAV. 

Contemporaneously  during  the  development  of  the  results  in  this  manuscript,  McGee, 
Spry  and  Hedrick11  have  addressed  the  problem  of  trajectory  planning  in  the  presence  of 
constant  wind.  The  mathematical  method  employed  by  them  appears  similar  to  ours  for  the 
case  of  constant  wind.  They  classify  the  set  of  initial  and  final  states  into  Types  1  and  2, 
where  Type  1  states  have  a  minimum  time  solution  while  Type  2  states  don’t.  Here  we  prove 
that  their  Type  2  states  form  a  set  of  measure  zero  which  means  that  they  are  sparse  in  the 
set  of  all  possible  states.  In  our  earlier  work,12  we  addressed  the  more  general  case  of  time 
and  space- varying  winds  and  proved  the  existence  of  time-optimal  solutions  (for  almost  every 
initial  and  final  states)  under  technical  conditions  involving  the  gradients  of  the  wind.  The 
verification  of  these  conditions  can  be  used  as  the  starting  point  for  a  numerical  algorithm  to 
compute  the  time-optimal  trajectory.  Here,  we  consider  time- varying  wind  vector  fields  (that 
is  W ( x ,  t)  —  W (t)  for  all  x  G  1R2  and  any  t  G  IR)  and  prove  the  convergence  of  an  algorithm 
for  computing  the  minimum-time  solutions.  Therefore,  the  work  in  this  manuscript  is  a 
generalization  of  the  work  of  McGee,  Spry  and  Hedrick.11 

Getting  back  to  the  equations  of  motion,  let  q  =  ( x The  optimal  control  problem  to 
be  solved  can  be  formulated  in  two  ways: 

Problem  1.  For  the  system  (1),  minimize  the  final  time  tj  subject  to  the  constraints: 
g(0)  =  qo  =  (iCtnVb)  and  q(tf )  =  q/  =  (xf^fd)',  «(•)  a  Lebesgue  measurable  function,  and 


satisfies, 


u(t)\  <  u 


max 


v 


Amin 


(2) 


for  all  t.  In  (2),  Rmin  is  the  minimum  turn  radius  in  the  absence  of  wind  and  arises  due  to 
mechanical  limitations  on  the  aircraft.  As  an  MAV  is  designed  to  fly  between  the  streets 
while  searching  for  targets,  we  must  ensure  that  the  MAV  enters  a  street  along  its  center 
line,  or  at  least  parallel  to  it.  In  our  numerical  examples,  iffd  is  the  angle  that  the  centerline 
makes  with  the  x-axis  of  the  coordinate  system.  Due  to  the  wind  vector  held,  the  actual 
tangent  line  to  the  solution  of  Problem  1  would  make  an  angle: 


Zi(f/)  =  arctan 


/  V  sin  ifi(tf)  +  W2(x,tf )  \ 
\V  cos  ifitf)  +  Wi(x,tf)  J 


fVsin  'ipfd  +  W2(x,tf)\ 
VH  cos  ipfd  +  W\{x,tf) )' 


(3) 


which  could  result  in  the  loss  of  the  MAV.  Hence,  it  is  clear  that  one  must  change  the 
end  condition  to  require  that  Zx(t/)  =  ipfd  to  minimize  the  possibility  of  aircraft  loss.  We 
denote  if(tf)  by  f°r  simplicity  and  also  denote  Wi(x,tf )  =  || PW (a; j, t/) || 2  cos Qw(xf,tf) 
and  W2(x,tf )  =  ||IV(a;/,f/)||2  sin 9w(xf,tf).  The  following  computations  yield  a  formula  for 
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V>/: 


V  sin-0/  +  W2(x,tf) 

V  cos  ip  (tf)  +  Wi(x,tf) 
V  sin  -0/  +  \\W(xf,tf)\\2  sin  0w(xf,tf) 
V cos 0/  +  ||VF(z/,£/)||2  cos 9w(xfltf) 

V  sin  (0/  -  ipfd) 

which  leads  to: 


=  tan  ipfd  (4) 

sin^fd  /5n 

cos  0/d 

=  \\W(xf,tf)\\2  sin(0/(i  -  9w(xf,tf)),  (6) 


0/  =  arcsin  sin(0/d  -  9w(xf,  t/))^  +  0/d;  (7) 


where 

0w(xf,  tf)  =  arctan  (8) 

This  leads  us  to  Problem  2: 

Problem  2.  For  the  system  (1),  minimize  tf  subject  to  the  constraints:  g(0)  =  qo  =  (xo,  0o); 
x(tf)  =  Xf\  and  Z x(tf )  =  0/^;  u(-)  a  Lebesgue  measurable  function,  and  satisfies  (2). 

In  this  paper,  we  will  focus  on  the  special  case  W (x,  t)  —  W ( t ).  An  existence  of  minimum¬ 
time  solution  result  for  the  more  general  case  of  time  and  space  varying  winds  (for  Problem 
1)  can  be  found  in  McNeely.12  In  the  same  work,  we  proved  the  uniqueness  of  minimum-time 
solutions  for  almost  every  combination  of  initial  and  final  states  for  Problem  1.  Here,  we 
prove  existence  and  uniqueness  of  minimum-time  solutions  for  the  case  of  time- varying  wind 
vector  fields  (that  is,  W(x,t )  =  W{t))  for  Problem  2. 

In  other  related  work,  Lissaman13  and  Zhao  and  Qi14  discuss  the  use  of  wind  gradients  to 
do  work  for  the  vehicle  while  in  flight.  Another  study  by  Venkataramanan  and  Dogan  is  based 
upon  the  use  of  a  nonlinear  controller  when  reconfiguring  the  formation  of  a  group  of  UAVs 
and  how  wind  effects  are  taken  into  consideration.15  Studies  by  Boyle  and  Chamitoff16  are 
based  on  a  stepwise  solution  to  a  constrained  optimization  problem  which  solves  a  nonlinear 
two-point  boundary  value  problem  at  each  timestep. 


The  Tour  Problem 

In  the  considered  mission,  a  small  unmanned  air  vehicle  (SUAV)  gives  the  MAV  a  set  of 
N  targets  X{  =  (xn ,Xi2);  i  =  1,  •  •  •  ,  IV  to  fly  over.  The  case  in  which  the  sensor  is  located 
on  the  aircraft  such  that  a  straight  line  flight  along  the  center-line  must  be  altered  in  order 
to  view  the  target,  is  outside  the  scope  of  this  paper  though  it  is  an  interesting  problem 
in  its  own  right.  For  each  target,  the  street  where  it  is  located  is  found  on  a  map  which 
yields  two  possible  angles  0,  or  0j  +  tt  for  flying  over  the  target.  The  mid-line  of  that 
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street  is  found  and  the  point  x*  is  projected  perpendicularly  from  the  target  onto  that  line. 
To  simplify  notation,  we  will  represent  this  new  point  also  as  X;  and  continue  to  work  with 
these  coordinates  henceforth.  For  the  i-th  target,  two  way-points  x\ ,  xf);  are  created  equally 
spaced  along  the  street  about  the  projected  point.  These  are  considered  entrance/exit  points 
for  viewing  the  target.  Thus,  each  exit  point  is  described  by  coordinates  q{  =  (x3il,xJi2,'i/jJi), 
where  i  —  1, 2,  ■  ■  •  ,  N  and  is  chosen  to  be  one  of  the  two  angles  ipi  or  +  n,  such  that  the 
unit  vector  (cos'ipj,  sin^)  points  away  from  the  target  (xi5  Ui)  at  the  point  Xj.  The  objective 
is  to  fly  from  target  i  to  target  k  ^  i  by  flying  from  a  point,  q\  to  qlk.  Thus  there  are  four 
possible  paths  corresponding  to  j,l  =  1,2.  In  this  paper,  we  do  not  address  the  problem  of 
flying  between  q\  to  q\  (that  is  flying  over  the  target  Xj)  where  r^s,  and  assume  it  to  be 
a  simple  straight  line  along  the  street.  Once  the  cost  of  flying  the  optimal  path  between  q\ 
to  qlk  for  all  i,  k  —  1,  2,  •  •  •  ,  TV;  i  ^  k  and  j,  l  —  1,2  then  one  can  set  up  a  task  assignment 
problem  as  outlined  in.3 

Discussion  of  Dubins’  Theorem  for  the  zero  wind  case 

We  briefly  discuss  Dubins’  result  on  the  existence  and  uniqueness  of  minimum  length 
solutions  for  a  particle  following  a  Lipschitz  continuous  path  at  a  constant  speed  V,  with  a 
bound  on  the  maximum  of  the  curvature  given  by  tt— •  This  theorem  states  that  for  every 
initial,  final  positions  and  velocities  the  minimum  time  (or  minimum  length)  solution  is  an 
arc-line-arc  or  arc-arc-arc  curve.  Here,  by  “arc”  we  mean  a  portion  of  a  circle  or  radius 
Rm in  and  “line”  means  a  straight  line  segment.  Furthermore,  for  the  arc-arc-arc  solution, 
Dubins’  theorem  states  that  the  middle  arc  must  have  length  greater  than  n  Rm-m .  For  the 
Tour  planning  problem,  Dubins’  theorem8  posits  the  existence  of  a  solution  to  the  minimum 
time  optimal  control  problem  for  the  special  case  W (x,  t)  =  (0,  0)  for  all  x  G  R2  and  t  6  1R+. 

As  the  minimum  time  solution  is  invariant  with  respect  to  translations  and  rotations  of 
the  coordinate  axis,  we  can  change  coordinates  so  that  the  initial  position  is  at  the  origin 
of  R2  and  the  final  position  is  along  the  x  axis.  The  initial  and  final  velocity  directions 
measured  with  respect  to  the  x-axis  are  termed  0o  and  0/  respectively  in  Figures  1-2. 
The  direction  of  the  velocity  induces  an  orientation  on  the  circles  tangent  to  the  velocity 
vector.  In  the  figures,  we  denote  the  center  of  the  counterclockwise  oriented  circle  by  Z$  and 
Zf  respectively,  while  the  centers  of  the  clockwise  oriented  curves  are  denoted  by  Yq  and 
Yf  respectively.  Thus  we  can  distinguish  between  the  Z0LY'f  from  the  Y'0LZf  arc-line-arc 
solutions  etc.  The  arc-line-arc  solution  is  found  when  the  initial  and  final  positions  are  fairly 
separated  as  in  Figure  1,  while  the  arc-arc- arc  solution  is  usually  found  when  the  initial  and 
final  positions  are  close  enough  as  in  Figure  2.  I  =  0  in  Figure  2  means  that  the  distance 
between  the  initial  and  final  positions  is  0.  Figure  2(a)  shows  a  Y0ZfYf  arc-line-arc  curve 
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connecting  the  initial  and  final  velocities.  Though  the  minimum  time  solutions  exist,  they 
need  not  be  unique  as  seen  in  Figures  1(b)  and  2(b).  A  more  important  point  from  the 
point  of  view  of  numerical  methods  is  the  differentiability  of  the  arc-lengths  as  a  function 
of  the  final  positions  and  velocities  while  holding  the  initial  positions  and  velocities  fixed. 
This  issue  is  discussed  in12  in  detail.  The  final  result  is  that  there  is  a  set  of  final  positions 
and  velocities  that  has  a  measure  zero,  and  where  the  minimum  time  curves  change  from 
arc-line-arc  to  arc-arc-arc  type  or  vice-versa  and  at  these  points  the  arc  length  function 
is  discontinuous  and  hence  non-differentiable.  It  is  also  possible  for  arc-line-arc  curves  to 
change  to  arc-line-arc  curves  in  a  discontinuous  manner. 

Time  optimal  trajectories  in  the  presence  of  wind 

For  a  non-zero  time-varying  (Caratheodory)  wind  vector  field,  it  is  possible  to  show 
existence  of  a  time-optimal  solution  if  we  make  the  following  assumptions.  Suppose  that  W 
is  a  differentiable  function  of  t  with  bounded  derivative  and  furthermore: 


Al.  l|W'|U  =  supt||W'(«)||<V'. 

A2.  ||fqu</Vr2-||tr||^«m„, 

where:  ||^y||oo  =  supt  ||^y(t)||2,  and  0  <  (3  <  1  is  some  constant. 

Remarks:  The  assumption  Al  is  easy  to  explain.  It  is  needed  to  ensure  that  the  MAV  will 
have  non- zero  forward  velocity  in  the  presence  of  the  worst-case  head  wind.  The  assumption 
A2  is  needed  because  the  MAV  has  a  limitation  on  the  maximum  rate  of  turn  given  by 
inequality  (2)  and  hence  cannot  keep  up  with  very  rapid  large  changes  in  the  wind  direction 
(please  refer  to  McNeely12  for  more  details). 

Solution  to  Problem  1 

In12  we  proved  the  following  theorem  (the  statement  has  been  modified  from12  to  include 
two  results  proved  in  the  same  manuscript): 

Theorem  1  Suppose  that  W(t )  is  a  time-varying  (Caratheodory)  vector  field  of  wind  ve¬ 
locities  satisfying  assumptions  Al  -  A2.  Then  there  exists  a  solution  to  Problem  1  for  the 
system  (1).  The  solution  is  unique  for  almost  every  collection  of  initial  and  final  states. 


Suppose  the  wind  W(x,t)  =  W(t)  depends  only  on  the  time  variable,  while  satisfying 
the  conditions  (Al)  -  (A2).  Then,  we  can  the  transform  coordinates  (with  the  same  time 
variable)  as  follows: 


x  =  <p(x,  t)  =  x  —  /  W(s)ds  and  = 


(9) 
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Let  q(t)  =  (a;(f),0(t))  be  the  solution  to  (1)  with  initial  condition  q0.  After  coordinate 
transformation,  this  solution  becomes: 

x(t)  =  x(t)  —  (  W(s)ds  and  -0(f)  = -0(f).  (10) 

Jo 

Therefore,  in  the  new  coordinates  {x,fi)),  the  system  equations  take  on  the  form: 

x  =  Id  (cos  0,  sin  0),  x(0)  =  x0]  (11) 

i  =  u,  -0(0)  =  0o,  (12) 


which  is  of  the  form  considered  by  Dubins.  The  final  states  of  the  transformed  system  for 
Problem  1  are  given  by: 

W(s)ds,  0(t/)  =  0/d.  (13) 

By  Theorem  1,  the  time  optimal  solution  exists  for  all  initial  and  final  conditions,  and  is 
unique  for  almost  every  set  of  initial  and  final  conditions  for  the  original  system.  Clearly,  the 
optimal  trajectory  must  be  a  Dubin’s  solution  for  the  transformed  system  with  initial  states 
(x,0)(O)  and  final  states  (x, 0)(t/)  (see  McNeely12  for  a  fuller  discussion).  The  problem  is 
that  tf  is  unknown  and  hence  x(tf )  is  also  unknown.  In  the  algorithm  in  Theorem  2,  we 
make  an  initial  guess  for  tf  and  improve  it  using  Newton’s  method.  The  convergence  of  the 
sequence  of  final  times  is  proved. 

Theorem  2  Select  >  0.  For  n  —  1,  2,  •  •  •  ,  consider  the  final  states: 

Tn~  1 

xir]-1)  =xf-  [  1  W(s)ds\  0(rj"1)  =  0/d.  (14) 

Jo 

Let  L(rj_1)  be  the  length  of  the  minimum-time  Dubins  trajectory  with  initial  state  (x,'0)(O) 
and  final  states  (x,  ^(rjC1),  and  denote 


x{tf )  =xf- 


T(r jT1) 


L_^n 

V 


(15) 


Set 


n  ^n—1 

Tf  =Tf 


.n—  1 
1 _ 

l  - 


(16) 


Then,  the  sequence  {rf]  n  G  7L}  converges  to  a  unique  rj  >  0  for  almost  every  set  of  initial 
and  final  conditions  for  the  original  system  provided  Tf  is  sufficiently  close  to  rf. 


9  of  33 


Proof:  Consider  the  function  /(r)  =  r  — T(r).  The  equation  /(r)  =  0  has  a  unique  solution 

Tf  by  Theorem  1  for  almost  every  set  of  initial  and  final  states.  This  is  because  /(rjf)  =  0 

yields  Tj  =  yr  .  This  means  that  the  minimum-time  trajectory  from  x(0)  to  the  point 
T*  T* 
xf  —  Jo  f  W(s)  ds  takes  exactly  time  rf  -  or,  in  other  words  -  x(rf)  =  Xf  —  f0f  W(s)  ds.  The 

existence  and  uniqueness  of  rf  having  been  addressed,  lets  consider  the  differentiability  of 

the  function  L(r).  As  studied  in  McNeely,12  this  function  is  continuously  differentiable  for 

almost  every  choice  of  r.  Hence,  provided  the  initial  selection  r9  is  sufficiently  close19  to  rf, 

Newton’s  method  is  applicable  to  the  solution  of  /(r)  =  0  which  results  in  (16). 

□ 

Remark:  Unfortunately,  it  is  difficult  to  characterize  how  “close”  the  initial  selection  r9  needs 
to  be  to  Tj  as  this  depends  on  xo,  Xf,  ifo,  Vh  and  the  wind  vectorfield  W{t)  in  a  complicated 
way. 

Solution  to  Problem  2.  We  will  first  prove  the  existence  of  time-optimal  solutions  for 
the  special  case  of  time- varying  wind  W(x,t)  =  W(t).  The  proof  for  the  general  case  also 
follows  the  same  pattern  though  the  equations  look  more  complicated.  As  the  algorithm  in 
this  paper  is  designed  for  the  case  for  W(x,t)  =  W(t),  this  is  the  significant  case  for  our 
purposes.  We  have  the  following  theorem: 

Theorem  3  Suppose  that  W(x,t )  =  W{t)  is  a  time-varying  (Caratheodory)  vector  field  of 
wind  velocities  satisfying  assumptions  Al  -  A2.  Then  there  exists  a  solution  to  Problem  2 
for  the  system  (1),  and  the  solution  is  unique  for  almost  every  collection  of  initial  and  final 
states. 

Proof: 

Existence:  The  proof  of  existence  proceeds  in  similar  fashion  to  the  proof  of  Theorem  1  as 
given  in.12  The  idea  is  to  establish  one  solution  and  then  use  Filippov’s  theorem12,18  to  prove 
the  existence  of  a  minimum-time  solution.  The  worst  case  wind  vector-field  is  the  one  that 
yields  the  system  equations: 

X  —  (V  +  ||VU||oo)  (cos-0  +  sin  VO;  =  u,  (17) 

as  the  instantaneous  turn-radius  given  by  R'min(t)  =  '  +lbUI°°  Rmin  is  the  largest 

possible.  The  condition  on  ip{tf)  becomes:  =  'iffd-  Dubin’s  Theorem8  is  applicable  to 

system  (17)  with  initial  states  (xojV’o)  and  final  states  and  it  yields  the  existence 

of  a  minimum-time  solution  for  the  worst-case  wind.  For  a  given  wind  field  W(t),  one  can 
construct  a  piecewise  continuous  u(t)  such  that  the  solution  stays  on  the  trajectory  for  the 
worst-case  wind  as  shown  in  the  proof  of  Theorem  1  in  McNeely.12  This  solution  satisfies 
Ax(tf)  =  ipfd  for  some  tf  >  0,  that  is  exactly  what  was  desired.  However,  this  solution  is 
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not  the  minimum  time  optimal  solution  for  the  wind  Wit),  and  we  prove  the  existence  of  the 
latter  using  Filippov’s  theorem.  In  order  to  apply  the  theorem,  we  need  to  make  a  change 
of  variables.  Let 

a(t)  =  arctan 

Then  we  obtain  (the  variable  t  is  suppressed  for  brevity): 


/  V sin  i/j(t)  +  W2(t)  \ 
cos +  W\(t)  J 


{V  cos'i/j  +  Wi)(V  ucos'ijj  +  W2)  —  (Id  sin -0  +  W2){— Vusmi/j  +  W\)  ,  . 

D2  +  ||ID||2  +  2D(IT2sin^  +  IFicosV0 

_  V2u  +  Vu  (W2  sin^  +  W1  cos^)  +  V  {W2  cos ^  -  W  sin^)  +  (W1  W2  -  W2  W) 

“  V2  +\\W\\2  +  2V  (W2  sin  ip  +  W\  cos  ijj)  (  ' 

The  above  equation  can  be  written  entirely  in  terms  of  a(t)  using  the  equations: 

ij)(t)  =  arcsin  (  ^  ^ )  1 1 2  _  (^(t))')  +  a(t);  Owif)  =  arctan  lyryT-  (21) 

V  L  )  Wi{t) 

The  system  in  variables  (x,  a)  has  initial  and  final  conditions  as  follows: 


x(0)  =  x0;  x(tf )  =  Xf 

V  sini/’o  +  W2{  0) 


ck(0)  =  arctan 


V  cos  '0o  +  VFi(O) 


a(tf)  =  ’hi 


(22) 

(23) 


It  is  straightforward  to  see  that  the  system  in  the  variables  (x,  a)  satisfies  all  the  conditions  of 
Filippov’s  theorem.12,18  As  we  have  already  shown  the  existence  of  one  solution  that  satisfies 
the  initial  and  final  conditions,  there  exists  a  measurable  control  u*(-)  and  a  minimum-time 
solution  x* (•),«*(•)  that  satisfies  the  boundary  conditions  for  almost  all  initial  and  final 
states. 

Uniqueness:  To  show  the  uniqueness  of  the  solution  for  almost  all  initial  and  final  states,  we 
need  to  work  with  the  x,^>  coordinates  as  in  (9): 

x  =  k^cos^,  sin-0),  x(0)  =  xo,  x(tf)  —  Xf—  /  W(s)ds;  (24) 

«,  =  u,  aretan  (vZlttj) t Wi) ) =  i’,d'  (25) 

By  our  proof  of  existence,  we  know  that  there  exists  a  minimum  time  solution  for  almost  all 
initial  and  final  states  for  the  original  system.  Assuming  it  exists,  let  (x,  -0)0),  t  G  [0,t/]  be 
this  solution.  Denoting  x(tf)  by  Xf  and  -0(tj)  by  -0^,  the  final  states  in  (24)  and  (25)  can  be 
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written  as: 


xf  =  xf  ~  [  W(s)ds ;  (26) 

Jo 

=  arcsin  sin (ipfd  -  6w(tf ))^  +  ipfd\  Ow(tf)  =  arctan  (27) 

With  (xf,ipf)  as  the  desired  final  state  for  the  transformed  system,  we  see  that  Dubins 
Theorem8  applies  and  there  exists  a  unique  solution  (x',ip,)(-)  defined  on  [0,  t*f]  for  almost 
every  combination  of  initial  and  final  states  of  the  transformed  system,  where  t*f  <  tf.  We 
can  transform  the  solution  back  to  the  original  coordinates  ( x ',ip')  while  remembering  that: 


x'(t})  =  Xf,  V't  *})  =  Pf, 


(28) 


to  obtain: 


=  x'(t*f)+  /  W(s)ds ; 
Jo 

f*} 

=  Xf  +  I  W(s)ds ; 

t} )  =  ipf- 


'tf 


(29) 

(30) 

(31) 


A  sufficient  condition  for  x'{t*f)  =  Xf  and  ip' iff)  =  ipf  is  t*j  =  tf.  If  t*f  <  tf,  then  it 
follows  that  both  solutions  ( x,ip)(t ),  t  e  [0,  t/]  and  (. x',ip')(t ),  t  e  [0,t^]  satisfy  the  boundary 
conditions,  and  the  first  solution  is  not  the  minimum-time  solution  contradicting  our  original 
assumption.  Hence  t*f  =  tf  is  both  necessary  and  sufficient  for  x'(t*f)  =  Xf  and  ip'(t*f )  =  ipf. 
As  the  Dubins  solution  in  the  transformed  coordinates  is  unique  for  almost  every  choice  of 
end  conditions,  the  same  is  true  for  the  minimum  time  solution  in  the  original  coordinates 
as  well. 

□ 

The  same  algorithm  and  convergence  proof  for  Problem  1  as  summarized  in  Theorem  2  also 
applies  to  Problem  2  with  a  minor  change. 
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Theorem  4  Select  r f  >  0.  For  n  =  1,  2,  •  •  •  ,  consider  the  final  states: 


—  (  n— 1\ 

[ 

x{Tf  ) 

=  xf-  / 

Jo 

v^r1) 

=  arcsin 

W(s )  ds; 


(32) 


y  sin  (fifd-6w(Tf  ))J+fifd',  where  (33) 

WFt?-1) 

(34) 


WfirT1) 

%(r;')=arctan— zjy. 


Lei  L(r)i  *)  6e  i/ie  length  of  the  minimum-time  Duhins  trajectory  with  initial  state  (x,fi)(0) 
and  final  states  (x,  '0)(rj”1),  and  denote 


T(r r1) 


L(r, 


n—  1  \ 


Lei 


_  _n-i  r/ 
r/  “  Tf 


V 

n—  1  rj-i/n— 1\ 


T(r 


/ 


1  - 


(35) 


(36) 


Then,  the  sequence  {rf]  n  G  converges  to  a  unique  Tj  >  0  for  almost  every  set  of  initial 
and  final  conditions  for  the  original  system  provided  t9  is  sufficiently  close  to  r*f. 


Proof:  The  proof  is  identical  to  that  of  Theorem  2,  except  for  the  fact  that  Theorem  3  is 
invoked  in  the  proof  instead  of  Theorem  1. 

□ 

Next,  we  describe  the  algorithms  given  in  Theorem  2  for  Problem  1,  and  Theorem  4  for 
Problem  2  in  greater  detail.  In  Step  2,  we  change  coordinates  for  numerical  stability  from 
(x,  if)  to  (z,  6)  such  that  in  the  new  coordinates  z0  is  at  the  origin  and  Zf  is  along  the  positive 
^i-axis.  One  could  also  employ  a  scaling  of  the  axes  in  this  step.  In  Step  3,  we  set  the  wind 
to  zero  and  use  Dubins  method  to  obtain  an  initial  estimate  for  the  final  time.  In  Step 
4,  we  modify  the  final  angle  6fd  to  Of  for  Problem  2.  This  step  is  not  necessary  for  Problem 
1.  At  Step  n,  the  estimate  to  tf  is  r”,  and  yields  Of.  On  convergence  of  the  algorithm  the 
correct  value  for  Of  will  be  found  along  with  tf. 

Algorithm  for  time  optimal  paths  in  the  presence  of  time-varying  wind 

1.  Start 

Input  q0  =  (xio,  x2o,  fi’o)  and  qf  =  (x\f,  x2/,  fifd)  {fi’fd  is  the  desired  final  angle),  mini¬ 
mum  turn  radius  ( Rm.in ),  speed  (' V ),  and  time-varying  wind  velocity  (W(t)).  Select  a 
tolerance  e  >  0. 
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2.  Perform  a  coordinate  transformation  for  numerical  stability 


Let  cos  0  =  nlf  Xl?,  and  sin 0  =  jf2/  X2?, 

r  If/— loll  |f/— ®o|| 


Define:  /?,  = 


cos  <p  —  sm  <p 
sin  0  cos  0 

We  use  the  rotation  matrix  i?  and  a  shift  of  origin  so  that  (x,  0)  is  transformed  into 
(z,  9),  with  Zf  along  the  positive  zi-axis.  Specifically,  z  =  RT  (x  —  xq),  and  0  =  0  —  0. 
We  have  to  consider  4  cases  to  compute  0  in  an  actual  implementation: 


(a)  If  cos  0  =  0  and  sin  0  >  0  let  0  =  7t/2 

(b)  If  cos  0  =  0  and  sin0  <  0  let  0  =  — 7t/2 

(c)  If  cos  0  >  0  let  0  =  arctan 

(d)  Else  0  =  arctan  +  tt 
Denote:  0O  =  0O  —  0  and  0/d  =  0/d  —  0. 

£o  and  Zf  are  computed  according  to  the  coordinate  transform:  z0  =  (0,0);  Zf  = 
Rt  (xf  —  Xo).  The  wind  vector  held  W(t)  is  transformed  to  W{t)  =  RT  W{t).  The  new 
system  of  equations  are: 


i  =  U(cos0,sin0)  +  W(t);  0  =  u.  (37) 

3.  Step  0:  Use  Dubins  method  to  find  an  initial  value  for 

Use  (z0,  0o )  and  ( Zf ,  0/)  to  compute  a  Dubins  trajectory  for  the  transformed  system 
with  zero  wind.  Let  L  be  the  length  of  this  trajectory.  Set  =  y. 

4.  Step  n:  Find  the  angle  9fd 

Find  0/  dependent  on  9fd,  V,  and  W{rf).  Figure  3  describes  the  idea  behind  the 
computation. 

Denote  by  V  =  (Vj,U2)  the  vector  D(cos0/,sin0/).  At  t  —  tJ,  Consider  V  + 
W(Tf)  =  (Vi  +  Wi(Tf),V2  +  W2(Tf))  where  V  is  as  yet  unknown  (because  rf  is  un¬ 
known)  . 

»T  i  ,,  ;  ,  a  v2+w2( rp 

Note  that  tan0/d  =  Vi+w^rp- 

Therefore 

Ui  +  W^r;)  =  \\V  +  W(t?)\\  cos0/d  (38) 

U2  +  W2{t7;)  =  \\V  +  W{t]) ||  sin0/d  (39) 
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If  we  combine  these  equations  and  solve  for  V2,  we  have  the  following  equation. 

V2  =  (V1  +  W1  (7?) )  tan  9fd  -  W2 (7?)  (40) 

We  also  have: 

Vl  +  V 22  =  \/2  (41) 

Solve  for  V\  and  V2  using  equations  (40)  and  (41). 

As  cos  (Of)  =  y  and  sin(61/)  =  y,  we  find  Of  using  the  following  process: 

(a)  If  cos  Of  =  0  and  sin  0 /  >  0  let  Of  =  7t/2 

(b)  If  cos  Of  =  0  and  sin  Of  <  0  let  Of  =  — 7t/2 

(c)  If  cos  Of  >  0  let  0t  =  arctan  sin  °J 

(d)  Else  Of  =  arctan  +  n 

V  /  J  cos  6f 

5.  Step  n  continued:  Compute  Tf 

T(n-1) 

Let  Zf  =  Zf  —  f0f  W(s)  ds.  Use  Dubins  method  with  initial  state  (zo,0o)  and  final 
state  ( Zf ,  0j_1).  Denote  the  length  of  the  path  divided  by  the  speed  V  by  T(rjn  ^)  as 
it  is  a  function  of  Tj" 

6.  If  | r-r}n"1)|  <  e,  then  go  to  Item  7.  Otherwise,  compute  the  solution  rj?-  to  /(r)  = 

r  (  n— 1\ 

T(t)  —  t  using  Newton’s  Method  or  Modified  Newton’s  Method.19  r?  =  rU1  —  r/n_i  . 

J  J  J  (Tyr  ) 

The  derivative  term  in  the  denominator  is  computed  numerically.  Go  to  Item  4. 

7.  Stop  We  take  Tf  as  the  approximate  value  of  Tf. 

Remark:  For  the  case  of  constant  wind  W(t)  =  W,  the  algorithm  is  particularly  simple  to 
implement  and  intuitive.  It  can  be  interpreted  as  follows.  We  seek  to  find  a  state  ( xj,^f ) 
such  that  in  the  absence  of  wind,  the  Dubins  solution  takes  time  tf  to  reach  this  state  from 
(xo,  V’o)-  The  property  of  this  state  is  that  it  takes  exactly  time  tf  for  a  particle  to  go  from 
Xj  to  Xf.  The  same  interpretation  also  holds  for  the  case  of  time- varying  wind. 

Numerical  Results 

To  illustrate  the  effectiveness  of  the  algorithm,  we  look  at  a  few  numerical  examples  while 
considering  a  constant  wind  vector  field.  I11  practice,  it  is  not  possible  to  have  knowledge 
of  the  wind  vector  field  W(t)  before  hand,  and  hence  constant  wind  is  the  important  case 
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to  consider.  We  first  consider  a  simple  trajectory  planning  problem  with  given  initial  and 
desired  final  positions  xo,  xj  and  headings  i/jq,  ipfd-  We  consider  Problem  2  for  several  cases 
with  various  constant  wind  vector  fields.  For  all  cases,  we  set  V  =  40  m/s  and  Rmin  =  100  m. 
The  units  of  distance  in  the  discussion  below  is  meters,  that  of  speed  is  meters  per  second, 
and  that  of  angle  is  radians. 

Trajectory  Planning  Example 

We  consider  two  generic  initial  and  final  states,  where  the  final  states  where  chosen  using 
the  random  number  generator  in  MATLAB.  The  wind  vectors  range  from  W  =  [—10;  —15] 
to  W  —  [10;  15]  and  all  satisfy  the  condition  \\W\\2  <  V.  Specifically,  the  states  were  q0  = 
(— 650,  — 100,  7t/4)  and  qf  =  (-692.9,433.62,-0.1824).  The  tolerance  e  for  the  algorithm 
was  chosen  to  be  e  =  0.1.  Applying  the  algorithm  for  Problem  2,  the  results  shown  in 
Figure  4  and  summarized  in  Table  1  were  obtained.  Figure  4  shows  the  actual  trajectories 
computed  and  their  variation  with  the  wind  can  be  seen.  Notice  that  all  trajectories  are 
tangential  to  each  other  at  the  end  point  in  spite  of  the  different  wind  conditions,  which  is 
exactly  as  required  for  solutions  to  Problem  2.  As  the  wind  increases  from  W  =  [10;  10]  m/s 
to  [10;  15]  m/s,  a  qualitative  change  in  the  trajectory  can  be  observed.  This  could  perhaps 
be  explained  by  observing  that  the  minimum  time  trajectory  would  be  the  one  with  the 
maximum  “coasting”  in  the  direction  of  the  wind. 

Next,  we  consider  a  case  where  the  initial  and  final  positions  are  identical  -  to  be  precise, 
qo  =  (0,0,0)  and  qf  =  (0,0 ,pi).  It  is  well-known  that  this  problem  has  a  non-unique  arc- 
arc-arc  solution  for  the  zero-wind  case  and  the  algorithm  computes  one  of  the  solutions  (the 
clockwise  one)  as  shown  in  Figure  5.  The  same  figure  also  shows  the  solution  to  Problem  2 
in  the  presence  of  constant  wind  W  =  [—5;  15]  m/s.  Notice  that  now  the  symmetry  is  lost, 
and  the  solution  to  Problem  2  is  a  clockwise  trajectory  which  is  tangential  to  the  x-axis,  as 
required. 

Optimal  Tour  Planning 

In  this  subsection,  we  consider  the  effect  of  the  wind  on  the  tour  planning  or  task  assign¬ 
ment  for  MAV’s.  Given  several  target  locations,  an  initial  starting  point  and  heading,  we 
would  like  to  find  the  order  in  which  to  visit  the  targets  such  that  the  total  time  of  travel 
is  minimized.  To  do  this,  we  first  compute  exit /entry  points  for  each  target  as  described  in 
the  section  on  the  Tour  Problem.  Then  we  set  up  a  directed  graph  containing  all  possible 
paths  between  the  targets  themselves,  and  between  the  starting  state  and  the  targets.  The 
minimum  time  taken  for  each  state  transition  is  computed  using  the  algorithm  and  assigned 
as  the  cost  of  the  state  transition.  Finally,  one  needs  to  search  the  cost  table  for  the  directed 
graph  to  find  a  state  transition  sequence  of  minimum  cost  that  starts  with  the  initial  state 


16  of  33 


and  covers  all  the  targets.  Several  publications1, 2,9  in  the  literature  have  addressed  the  search 
problem  and  that  is  not  the  topic  under  investigation  in  this  paper.  As  our  interest  is  in 
demonstrating  the  effect  of  wind  on  the  task  assignment  problem,  we  consider  two  examples 
with  3  targets  each  and  the  direct  graph  search  was  done  manually  with  some  effort.  The 
wind  for  both  examples  was  set  at  [10;  15]  m/s  though  this  is  just  a  number  and  the  theory 
developed  earlier  does  not  depend  on  any  specific  wind  speeds. 

In  the  first  example,  the  drop  off  point  is  (—650,  —100,  |),  and  in  the  second  example, 
the  drop  off  point  is  (400,400,  ^),  and  they  have  been  denoted  by  S  in  Figures  6(a)  and 
6(b).  The  three  targets  are  Target  1  =  (—703,434.9),  Target  2  =  (790.9,632.2),  and  Target 
3  =  (-179.5,-124.8).  Note  that  these  targets  were  chosen  using  the  map  of  a  simulated 
city  that  is  available  in  the  MultiUAV  simulation.3 

The  entry/exit  points  for  each  of  the  targets  were  found  as  discussed  in  Section  on  the 
Tour  Problem.  Those  about  the  first  target  are  la  =  (-682.9,433.6,-0.18)  and  16  = 
(—722.2,440.9,2.96).  The  second  target  has  bounds  2 a  =  (810.8,629.8,-0.19)  and  26  = 
(771.53,637.19,2.9565).  The  third  target  has  entry /exit  points  3a  =  (—158.5,-121.0,-0.18) 
and  3c  =  (-197.8,-113.7,2.96).  For  the  first  example,  the  Tables  2  and  3  were  obtained 
for  the  state  transition  times,  where  the  time  is  measured  in  seconds.  The  time  optimal 
tours  are  plotted  in  Figure  6(a).  The  results  of  the  minimum-time  tour  computation  showed 
an  improvement  in  the  time  needed  to  travel  the  minimal  path  in  these  wind  conditions. 
The  optimal  path  where  there  is  no  wind  present  was  computed  to  be  S-lb-la-3b-3a-2b-2a  - 
which  is  travel  from  the  drop  off  point  to  the  first  target,  then  the  third  target,  and  finally  to 
the  second  target.  The  cost  of  traveling  this  path  in  the  case  of  no  wind  was  69.4  seconds.  In 
the  presence  of  wind,  the  optimal  path  was  found  to  be  S-la-lb-3b-3a-2b-2a  and  the  cost  got 
reduced  to  64.9  seconds.  Although  travel  to  each  target  was  performed  in  the  same  order, 
the  approach  was  different  in  each  scenario.  This  case  shows  how  the  presence  of  wind  could 
help  the  trajectory  of  an  MAV. 

For  the  second  example,  only  the  relationships  with  S  and  the  other  points  change.  The 
cost  tables  for  this  problem  are  given  in  Tables  4  and  5.  The  time  optimal  tours  for  the  zero 
wind  and  wind  case  is  shown  in  Figure  6(b).  The  results  of  the  search  routine  showed  that 
the  minimal  path  in  wind  conditions  was  longer  than  that  in  no  wind.  The  optimal  path 
where  there  is  no  wind  was  S-2a-2b-3a-3b-la-lb  with  cost  75  seconds.  In  the  presence  of 
wind,  the  optimal  path  was  S-3a-3b-la-lb-2b-2a  and  the  cost  was  92.8  seconds.  This  path 
is  very  different  from  that  where  no  wind  is  present. 
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Conclusion 

In  this  paper,  we  consider  the  problem  of  constructing  minimum  time  trajectories  for  a 
Dubins  vehicle  in  the  presence  of  a  time  varying  wind  vectorfield.  We  have  presented  results 
on  the  existence  and  uniqueness  of  minimum-time  solutions  for  a  Dubins  vehicle  flying  in  a 
general  time-varying  wind  vectorfield  under  some  technical  conditions.  Using  these  results 
we  have  proposed  and  proved  the  convergence  of  an  algorithm  for  obtaining  minimum-time 
solution  for  such  a  vehicle.  We  have  presented  the  results  of  numerical  experiments  that  show 
that  the  importance  of  considering  wind  vectorfields  while  planning  the  tour  for  a  Dubins 
vehicle. 
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Table  1: 
Table  2: 
Table  3: 
Table  4: 
Table  5: 


List  of  Table  Captions 

Wind  vectors  and  minimum  times  for  optimal  trajectories. 

Tour  example  1:  zero  wind  case. 

Tour  example  2:  wind  case. 

Tour  example  2:  zero  wind  case. 

Tour  example  2:  wind  case. 
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Wind  vector 
(m/s) 

Time 

(s) 

[-10,-15] 

23.33 

[-5,-10] 

19.75 

[0,0] 

16.61 

[5,5] 

17.22 

[10,10] 

20.25 

[10,15] 

19.11 

Table  1: 


22  of  33 


S 

la 

lb 

2a 

2b 

3a 

3b 

s 

- 

14.6 

17.3 

45.8 

40.2 

21.4 

11.6 

la 

- 

- 

1 

44.3 

36.7 

24.5 

18.6 

lb 

- 

1 

- 

52.6 

44.3 

32.2 

24.5 

2a 

- 

44.3 

52.6 

- 

1 

35.0 

40.6 

2b 

- 

36.7 

44.3 

1 

- 

30.6 

35.0 

3a 

- 

24.5 

32.2 

35.0 

30.6 

- 

1 

3b 

- 

18.6 

24.5 

40.6 

35.0 

1 

- 

Table  2: 
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S 

la 

lb 

2a 

2b 

3a 

3b 

s 

- 

11.5 

19.6 

35.9 

29.5 

15.6 

10.6 

la 

- 

- 

0.9 

35.4 

29.6 

26.8 

24.4 

lb 

- 

1.2 

- 

43.0 

37.2 

31.5 

28.3 

2a 

- 

66.3 

72.7 

- 

0.9 

60.7 

65.5 

2b 

- 

57.4 

63.6 

1.2 

- 

54.4 

58.5 

3a 

- 

29.0 

37.0 

27.0 

21.8 

- 

0.9 

3b 

- 

18.1 

26.6 

32.1 

25.6 

1.2 

- 
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S 

la 

lb 

2a 

2b 

3a 

3b 

s 

- 

27.3 

36.1 

22.8 

19.8 

19.5 

23.9 

la 

- 

- 

1 

44.3 

36.7 

24.5 

18.6 

lb 

- 

1 

- 

52.6 

44.3 

32.2 

24.5 

2a 

- 

44.3 

52.6 

- 

1 

35.0 

40.6 

2b 

- 

36.7 

44.3 

1 

- 

30.6 

35.0 

3a 

- 

24.5 

32.2 

35.0 

30.6 

- 

1 

3b 

- 

18.6 

24.5 

40.6 

35.0 

1 

- 
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S  la 

lb 

2a 

2b 

3a 

3b 

s 

-  39.6 

46.4 

19.1 

18.3 

34.5 

38.6 

la 

- 

0.9 

35.4 

29.6 

26.8 

24.4 

lb 

-  1.2 

- 

43.0 

37.2 

31.5 

28.3 

2a 

-  66.3 

72.7 

- 

0.9 

60.7 

65.5 

2b 

-  57.4 

63.6 

1.2 

- 

54.4 

58.5 

3a 

-  29.0 

37.0 

27.0 

21.8 

- 

0.9 

3b 

-  18.1 

26.6 

32.1 

25.6 

1.2 

- 
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List  of  Figure  Captions 

Figure  1:  Illustration  of  arc-line-arc  solutions  for  the  minimum  time  optimal 
control  problem. 

a) :  Illustration  of  an  arc-line-arc  Y^LZf  minimum  length/time  solution  for 
sufficiently  separated  initial  and  final  positions. 

b) :  Illustration  of  non-unique  arc-line-arc  minimum  length  solutions. 

Figure  2:  Illustration  of  arc-arc-arc  solutions  for  the  minimum  time  optimal 
control  problem. 

a) :  Illustration  of  an  arc-arc-arc  solution  when  the  initial  and  final  positions 
coincide. 

b) :  Illustration  of  non-unique  arc-arc-arc  solutions. 


Figure  3: 

Relationship  of  Ofd  and  Of. 

Figure  4: 
ditions. 

Minimum  time  arc-line-arc  solutions  for  zero  and  non-zero  wind  con- 

Figure  5: 
ditions. 

Minimum  time  arc-arc-arc  solutions  for  zero  and  non-zero  wind  con- 

Figure  6: 
targets. 

Minimum  path  examples  comparing  wind  and  no  wind  tours  of  three 

a) :  Tour  planning  problem  1. 

b) :  Tour  planning  problem  2. 
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Illustration  of  an  arc-line-arc  Y0LZf  minimum  length/time  solution  for  sufficiently 

separated  initial  and  final  positions. 


Illustration  of  non-unique  arc-line-arc  minimum  length  solutions. 


Figure  1: 
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Illustration  of  an  arc-arc-arc  solution  when  the  initial  and  final  positions  coincide. 


Illustration  of  non-unique  arc-arc-arc  solutions. 


Figure  2: 
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V(cos  Of,  sin  Of) 


Figure  3: 
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Figure  4: 
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Figure  5: 
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Figure  6: 
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